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Abstract 

The effect of neutrino trapping on the longitudinal dielectric function at low densities has been 
investigated by using different relativistic mean field models. Parameter sets G2 of Furnstahl- 
Serot-Tang and Z271 of Horowitz-Piekarewicz, along with the adjusted parameter sets of both 
models, have been used in this study. The role of the isovector adjustment and the effect of the 
Coulomb interaction have been also studied. The effect of the isovector adjustment is found to 
be more significant in the Horowitz-Piekarewicz model, not only in the neutrinoless matter, but 
also in the matter with neutrino trapping. Although almost independent to the variation of the 
leptonic fraction, the instability region of matter with neutrino trapping is found to be larger. The 
presence of more protons and electrons compared to the neutrinoless case is the reason behind 
this finding. For parameter sets with soft equation of states at low density, the appearance of 
a large and negative £l(q,Qo = 0) in some parts of the edge of the instability region in matter 
with neutrino trapping is understood as a consequence of the fact that the Coulomb interaction 
produced by electrons and protons interaction is larger than the repulsive isovector interaction 
created by the asymmetry between the proton and neutron numbers. 

PACS numbers: 13.15.+g, 25.30.Pt, 97.60. Jd 



1 



I. INTRODUCTION 



At low densities both the relativistic and the non-relativistic mean field models predict a 
liquid-gas phase transition region for nuclear matter leading, for dense star matter, to a non- 
homogeneous phase commonly called pasta phase, which is formed by a competition between 
the long range Coulomb repulsion and the short range nuclear attraction l| . This transition 
has substantial consequences on the properties of stellar matter and neutrino transport 
Considerable efforts to comprehend the uniform ground state stability of multi-component 
systems consisting of electrons, neutrinos, protons, and neutrons as a good approximation of 
this transition have been recently devoted, not only in the zero temperature approximation, 



but also for finite temperature 



y devoted, not only n 
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12 1 . It is obvious that in order to 



understand the physics inside the non-homogeneous (unstable) regions like the mechanism 
of nuclear creation with slab-like or rod-like shape, we have to go beyo nd the mean field 
approximation. Attempts in this direction are discussed in Refs. Moreover, 
in the collapsing supernova core and at sub-nuclear densities, the transition of nuclear shape 
from sphere to other exotic shapes has significant effects to the neutrino mean free path. 
However, how these effects modify the neutrino mean free path is not fully understood 



yet 
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Another motivation to study this transition comes from the fact that a neutron star 
is expected to have a solid inner crust of nonuniform neutron-rich matter above its liquid 

n 

mantle [5j and the mass of its crust depends sensitively on the density of its inner edge and 
on its equation of state (EOS) [4j]. On the other hand, the critical density (p c ), a density 
at which the uniform liquid becomes unstable to a small density fluctuation, can be used as 
a good approximation of the edge density of the crust 5|. By generalizing the dynamical 
stability analysis of Ref. 17] in order to accommodate the various nonlinear terms in the 
relativistic mean field (RMF) model of Horowitz-Piekarewicz jg], Carriere et al. [5] found 
a strong correlation between p c in the neutron star and the density dependence of nuclear 
matter symmetry energy (a sym ). This leads to a suggestion that a measurement of the 
neutron radius in 208 Pb will provide useful information on the p c 

In our previous work 11], the critical densities of uniform matter with and without 
neutrino trapping have been calculated and analyzed by means of different RMF models. In 
this analysis it is shown that the interplay between the dominant contribution of the matter 
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composition and the effective masses of mesons and nucleons leads to higher critical densities 
for matter with neutrino trapping. Furthermore, it was also found that the predicted critical 
density is insensitive to the number of trapped neutrinos as well as to the RMF model used. 
However, the discussion about the reason behind these findings was not quite robust. On the 
other hand, as we mentioned above, the neutrino transport is very crucial in the dynamics 
of the core-collapsing supernova due to the fact that the neutrinos carry most of the energy 
away and will lose also their energies by exciting collective nuclear and plasmon modes [8(. 
Similar situation can be also found in the neutron stars. Moreover, it was also shown in 
Ref. 8] that the behavior of electrons in matter depends strongly on the wave length or 
momentum of the external perturbation q. Note that this momentum is related to the 
energy transfer of the neutrinos that propagate in matter. 

The present paper reports on the extension of our previous investigation by calculat- 
ing the longitudinal dielectric function of ERMF models and analyzing the relation between 
the obtained results and the isovector sector adjustment, the presence of the long-range 
Coulomb interaction, as well as the presence of electrons. The purpose of this work is to 
explain the reason behind the appearance of each point along the onset of the instability. 
To this end, we should emphasize here that we need to calculate the dielectric function in 
the edge of non-homogeneous regions because information from the critical density alone is 
insufficient. Furthermore, it is also important to emphasize that our definition of the in- 
stability is not the non-homogeneous area, but rather it is connected with the points where 
these non-homogeneities start to appear. As a consequence, the assumption of the uniform 
matter in the calculation is still valid. 

This paper is organized as follows. The RMF models and some constraints used in the 
present analysis are briefly discussed in Sec. [Til In Sec. PH a discussion of the longitudinal 
dielectric function is given. In Sec. IIVI we present the graphical results of the onset of 
instability along with the corresponding discussions. Finally, we give the conclusion in 
Sec. ED 



II. RMF MODELS 



To describe the multi-component matter, we use the Lagrangian density 11] 

C = Cn + Cm + £hp + £l , 
3 



(1) 



where the first three terms describe the nucleons in the mean field level, while the last 
term indicates a free Lagrangian for leptons. The first term is the Lagrangian for nucleons 
interacting with each other via meson exchanges. The second term is the Lagrangian for 
mesons, containing also their nonlinear self coupling information. The third term is added 
to accommodate the Horowitz-Piekarewicz isovector nonlinear term [6|]. 

In this study we investigate two RMF models, namely, the G2 parameter set of Furnstahl- 
Serot-Tang 19( (also known as the ERMF model) and the Z271 parameter set of Horowitz- 
Piekarewicz |6j. The effective Lagrangian density of Furnstahl-Serot-Tang model has been 
constructed to fulfill the symmetries of quantum chromo dynamics and is expanded in powers 
of the fields and their derivatives up to order v = 4. Furthermore, the Lagrangian exploits 
the natural coupling constants and its application to study the properties of finite nuclei 
is quite successful. The inclusion of the higher order terms was found to be unimportant 
(and undetermined) for the nuclear observables of interest [191 ] . On the other hand, the 
Horowitz-Piekarewicz model is an extension of the standard RMF model with an additional 
isovector-vector nonlinear term. 

Therefore, the two models can be considered as the generalization of the standard RMF 
models. Details of the individual terms and coupling constants of both models can be found 
in Ref. ll|. For each model we use two different parameter sets, i.e., the standard ones (G2 



and Z271) and the adjusted ones which produce softer symmetry energy predictions at high 
densities (G2* and Z271*). Parameter set Z271* is obtained by adding an isovector-vector 
nonlinear term with the coupling constant Ay in the Lagrangian density of Z271 parameter 
set and followed by an adjustment of the g p and Ay [6]. On the other hand the G2* 
parameter set is obtained by using a similar procedure. However, since the ERMF model 
already contains an isovector-vector nonlinear term, the G2* parameter set is obtained by 
merely an adjustment of the g p and rj p parameters, keeping the symmetry energy at the 
same value with the G2 parameter set, i.e., E sym = 24.1 MeV at kp = 1.14 fm _1 11]. 

By comparing the low density instability regions for matter with and without neutrino 
trapping obtained from both parameter sets we can investigate the role of isovector terms 
and the correlation between the instability regions and the symmetry energy. Furthermore, 
in this study we use the zero temperature approximation. We note that in the real situation 
the temperature of protoneutron stars is not equal to zero and a supernovae inner core can 
have a temperature around T ~ (10-50) MeV. Indeed, the stability of uniform matter is 
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sensitive to temperature. This indicates that investigations at finite temperature will need 
to be addressed in the future. 

In our approximation, the following constraints can be used to determine the fraction of 
every constituent in matter which are later used to calculate the Fermi momentum of every 
constituent involved: 

• the balance equation for the chemical potentials 

tin + Pv e = Pp + Pe , (2) 

• conservation of the charge neutrality 

Pe = P P , (3) 

• and fixed electronic-leptonic fraction 

Y le =Y e + Y Ve , (4) 

where the total baryon density is limited by 

Pb = Pn + P P ■ (5) 

Note that in the case of matter without neutrino trapping we have Y Ve = and the value of 
Y\ is not fixed. 

I 

As has been reported in the previous work [ll||, besides the EOS, at low density regimes 
the asymmetry between the proton and neutron number (a=Y^ — Yp) in matter with and 
without neutrino trapping (NT) is also different, i.e., a wlthont NT is closer to the asymmetry of 
the pure neutron matter (PNM), whereas a wlth NT is closer to the symmetric nuclear matter 
(SNM). Thus, a without NT has a strong correlation with a sym , in contrast to the a with NT . 
This behavior is also found in the properties of Fermi momentum of each constituent in 
matter. Another different phenomenon is that in matter with NT we have k e F = k F ~ k F , 
while matter without NT has k e F = k F -C k F . The latter indicates that the role of isovector 
contribution is more significant in matter without NT than in matter with NT, while the 
Coulomb interaction has a more significant effect in matter with NT due to the presence 
of more protons and electrons. Note that Fermi momentum of every constituent is one of 
the required information, besides the nucleon effective mass, to calculate the polarizations 
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FIG. 1: (Color online) Proton (electron) Fermi momentum as a function of the ratio between 
baryon and nuclear saturation densities for the G2, G2*, Z271, and Z271* parameter sets in the 
neutrinoless matter and in matter with neutrino trapping with Yi e = 0.3 and 0.4. 

in the longitudinal dielectric function. In Fig. [TJ the effects of the neutrinos in matter on 
the electron or proton Fermi momentum at low density regimes are shown. 

The symmetry energies (a sy m) of the corresponding parameter sets are shown in Fig. [2j It 
is clearly seen that, different from the Horowitz-Piekarewicz model, the high density adjust- 
ment in the isovector- vector channel of the Furnstahl-Serot-Tang model does not significantly 
affect its a sym at low density. On the other hand, parameter set with softer symmetry energy 
at high density of Horowitz-Piekarewicz model becomes stiffer at low density regimes. It 
means that, in contrast to the Furnstahl-Serot-Tang model, the high density isovector ad- 
justment in the Horowitz-Piekarewicz model leads to a more repulsive isovector interaction 
than the standard one at low density regimes. The different forms of the nonlinear terms in 
both models are responsible for these different behaviors. 

From both figures, it can be explicitly seen that the presence of neutrinos in matter 
leads to higher Fermi momenta of protons and electrons compared to the case of matter 
without NT, while in the latter Fermi momentum of every constituent, which is represented 
by protons and electrons, is independent to the model used. Furthermore, in matter without 
neutrino trapping, due to its proton-neutron asymmetry that is closer to PNM, the Fermi 
momentum of each constituent can be correlated to the a svm - 
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FIG. 2: (Color online) Symmetry energy as a function of the ratio between baryon and nuclear 
saturation densities for the G2, G2*, Z271, and Z271* parameter sets. 

III. LONGITUDINAL DIELECTRIC FUNCTION 

The longitudinal dielectric function can be written as jlj] 



The uniform ground state system becomes unstable to small-amplitude density fluctuations 
with perturbation momentum q when EL(q,q = 0) = El < 0. Note that in Eq. (jSJ) qo is the 
time-component of the four-momentum q^ = (qo,q) and q — \q\. The critical density p c is 
the largest density for which the above condition has a solution. For matter consisting of 
protons, neutrons, and electrons, the longitudinal meson propagator is given by 



where d svp — —(d sv + d sp ), d svp — —{d sv — d sp ), d vp — d v — d p , g? 33 — d g + d v + d p + 2d vp and 
c?44 = d v + d p — 2d vp . In this form, mixing propagators between isoscalar-scalar and isoscalar- 
vector (d sv ), isoscalar- vector and isovector- vector (d vp ), isoscalar-scalar and isovector- vector 
{d sp ) are present due to the nonlinear mixing terms in the model, in addition to the standard 
photon, omega, sigma and rho propagators (d g , d v , d s and d p ). These propagators are 



£l(<7, qo) = det [1 - D L (q, qo)U L (q, q )\ . 



(6) 



f d g -d g \ 

— d s d~^ vp d svp 

~d g d svp d 33 d~ p 

\ d ivp d vp rf 44 / 



(7) 



determined from the quadratic fluctuations around the static solutions which are generated 
by the second derivatives of energy density (d 2 e/dfad(f)j), where fa and fa are the involved 
meson fields. The explicit forms of the a, u, and p propagators are 

, 9l(q 2 + m*J)(g 2 + m; 2 ) 

(q 2 + m* 2 )(q 2 + m* 2 )(q 2 + m* 2 ) + (II°J 2 (g 2 + m* 2 ) + (II° p ) 2 (g 2 + m* 2 ) ' U 

, 9l{q 2 + rnl 2 ){q 2 + m* p 2 ) 

(q 2 + ml 2 )(q 2 + m* p 2 )(q 2 + m* 2 ) + (II°J V + m£ 2 ) - + m* 2 ) ' U 
l/4^(g 2 + m; 2 )(g 2 + m: 2 ) 

p (g 2 + ml 2 )(q 2 + m* p 2 )(q 2 + m* 2 ) + (Tl» p ) 2 (q 2 + m* 2 ) - (II°°)V + m* 2 ) ' 1 j 

and in the mixing propagators, 

^n° g ( 9 2 + m; 2 ) 
#(g,g = o) ' 

l/2^X p (g 2 + < 2 ) no , 
H(q,q = 0) ' ^ 

l/2g p ^n°° p (g 2 + ml 2 ) 

^ = = ' (13) 

with 

H(q,q = 0) = (q 2 + m: 2 )(q 2 + m; 2 )(q 2 + m: 2 ) + (Ul w ) 2 (q 2 + m; 2 ) 

+ (nyV + < 2 )-(nS)V + < 2 ), (14) 

where the effective mass of each meson is 

< 2 = = m l + 2b2(T + 3b ^ 2 - d ^ V o - ^sb 2 , (15) 
ml 2 = - J^- = m 2 + 2rf 2 a + rf 3 a 2 + 3c 3 l/ 2 + A„6 2 , (16) 

™; 2 = = m l + 2 ^ + A ^ 2 + A ^ 2 , (17) 



and the effective mixing masses read 



nl = -Q^y = 2d 2 V + 2d 3 aV , (18) 
n °" = "9^ = 2/A + 2A * <r '' t, • (19) 



whereas the propagator of photon (Coulomb) is 



d 9 = J 2 - (21) 
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The longitudinal polarization matrix given in Eq. ([6]) reads 

/ n e 00 o o o \ 



n 



L 



o n s w m n» 
o ng o 
I o n: o n;/ 



(22) 



The formulas for polarization elements in are given in, e.g., Ref. [111 ]. Note that for 
the Horowitz-Piekarewicz model II° W and n° equal to zero. On the other hand, for the 
Furnstahl-Serot-Tang model, 11° and A s in n° p are zero. In Figs. [3] and HI we show some 
plots as examples of the polarizations and propagators at q = 0.5 fm _1 . It is clearly seen 
that all polarizations depend strongly on whether the neutrinos are trapped or not in matter. 
In matter with NT, at a density close to 0.6 pB, they have only a weak dependence on the 
model used due to the increasing role of the nucleon effective mass. For matter without NT, 
a similar model dependence on its Fermi momentum fc^ shows up (see Fig. [p. This indicates 
that the polarizations can be correlated with q syTri . Note that the nucleon effective mass has 



no correlation with the symmetry energy a sym [ll|. This demonstrates that, instead of the 
effective mass of nucleons, Fermi momentum of each constituent controls the behavior of 
each polarization. 

Figure H] shows the behavior of the propagators as a function of the ratio between baryon 
and nuclear saturation densities. The scalar (d s ), vector (d v ) and rho (d p ) propagators are 
clearly model dependent and they are insensitive to the presence of the neutrino (NT) in 
matter. Only d p depends on the isovector adjustment. The presence of neutrinos has a 
significant effect only in the mixing propagators (d sp and d vp ). However, their contributions 
are smaller compared to the rest. We also note that the propagator d sv , which appears in 
the Furnstahl-Serot-Tang model, has an order of magnitude that equals to the propagators 
d sp and d vp , and, furthermore, it does not depend on the adjustment of the isovector sector. 
These properties are a manifestation of the interplay among a, uj and p meson effective 
masses as well as mixing effective masses 11° 11° and 11° 

To simplify our discussion, let us neglect the minor contributions of the mixing prop- 
agators. In this picture we obtain d mth NT ~ rf wlthout NT . Together with the fact that d p 
depends on the adjustment of isovector sector (model dependent), this indicates that d p can 
be related with a sym . 
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FIG. 3: (Color online) The proton-scalar, -longitudinal, -mixing, and electron- longitudinal polar- 
izations as a function of the ratio between baryon and nuclear saturation densities for the G2, G2*, 
Z271, and Z271* parameter sets in neutrinoless matter and matter with neutrino trapping with 
Y\ e = 0.4. All curves have been obtained by using q = 0.5 fm _1 . 

Therefore, in the case of matter with NT, we can understand that the critical density is 
higher and the instability region is larger than those in the case without NT. Furthermore, 
the fact that both of them are insensitive to the value of Yi e and they are not influenced by 
the isovector treatment appears as a consequence of the small proton-neutron asymmetry 
a, which makes the role of the polarizations more dominant compared to the role of the 
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FIG. 4: (Color online) The scalar, vector, rho and mixing scalar-rho or vector-rho propagators as 
a function of the ratio between baryon and nuclear saturation densities for the G2, G2*, Z271, and 
Z271* parameter sets in neutrinoless matter and matter with neutrino trapping with Y[ e = 0.4 and 
q = 0.5 fm -1 . Note that for the Horowitz-Piekarewicz model d sp = 0, whereas for the Furnstahl- 
Serot-Tang model d vp =0. 

corresponding propagators, in controlling the behavior of the longitudinal dielectric function. 
The correlation between critical density, as well as the onset of the instability, and a sym in 
matter without NT are caused by two sources: the dependence of the polarization on the 
a sym as well as on the isovector propagators (d p , d sp and d vp ). 
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FIG. 5: (Color online) Onset of the instability of the neutrinoless matter obtained by using the 
Horowitz-Piekarewicz (Z271 and Z271*) and Furnstahl-Serot-Tang (G2 and G2*) models as func- 
tions of the ratio between baryon and nuclear saturation densities and the perturbation momen- 
tum q. 

IV. THE ONSET OF INSTABILITY 

As mentioned in the Introduction, here we intend to investigate every point in the onset of 
the instability region in a great detail. To this end, we plot the projection of the longitudinal 
dielectric function given in Eq. (jHJ) on the p/po — q plane in the case that El = 0. To obtain 
more information on the role of the Coulomb interaction in the limit of q ~ (almost zero 
perturbation), we also present the longitudinal dielectric function at q close to zero, i.e., 
£ h(q = 0.01 fm _1 ,g = 0), as a function of the p/po- 

The effects of the isovector-vector channel adjustment in both models are exhibited in 
Fig. Obviously, from the size and the position of the boundary of the instability region, 
the adjustment in the isovector-vector sector has a more significant effect in the Horowitz- 
Piekarewicz model compared with the Furnstahl-Serot-Tang one. It is also clear from the 
figure that this adjustment leads to a higher critical density in both models. This result is 



certainly consistent with our previous study (cf. Fig. 8 of Ref. |ll|] in the case of Y\ £ = Y e ). 
Another important finding obtained from these figures is that in both cases the ERMF model 
of Furnstahl-Serot-Tang yields a smaller onset of instability. This is due to the fact that the 
parameter sets Z271* and G2* have larger symmetry energies compared to the Z271 and G2. 
Therefore, beta-equilibrium of these parameter sets is attained with larger proton fractions 
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FIG. 6: (Color online) Onset of the instability of the neutrinoless matter obtained by using the 
Furnstahl-Serot-Tang model with and without (symbolized with — ) Coulomb interaction and elec- 
trons contributions. Note that the solid (red) and dash-dotted (blue) curves [as well as the dashed 
(green) and dotted (black) ones] are coincident. 
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FIG. 7: (Color online) Longitudinal dielectric function for matter without NT obtained by using 
the Furnstahl-Serot-Tang model with and without (symbolized with — ) Coulomb interaction as well 
as electrons contributions. Note that the dashed (green) and dotted (black) curves are coincident. 
All curves have been calculated by using q = 0.01 fm _1 . 

and a smaller contribution from the repulsive isovector channel. Thus, as we expected, Fig. [5] 
shows explicitly the correlation between the onset of the instability region with a sym . The 
reason of this fact has been explained in the previous section. 
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FIG. 8: (Color online) Onset of the instability for matter with NT as functions of the ratio between 
baryon and nuclear saturation densities and the perturbation momentum q. All curves have been 
obtained by using Yj e = 0.3. The results are obtained by using the Horowitz-Piekarewicz (Z271 
and Z271*) and Furnstahl-Serot-Tang (G2 and G2*) models. 

In Fig. [HI we show the effect of the electron absence (indicated by "— e" in the figure) on 
the onset of instability by switching off their contribution in the case of matter without NT 
using G2 parameter set, with and without Coulomb contribution. Note that the latter is 
indicated with "—Coulomb" in the figure. The effect of the electron presence on the size of 
the instability region is found to be negligible in both cases. The reason is that the number 
of the electrons is too small in matter without NT, and therefore, the effect of the attractive 
Coulomb interaction generated by electrons and protons is too weak to produce a visible 
impact on every point in the onset of instability. Such behavior is observed even in the limit 
of q close to zero (see Fig. [7]). 

On the other hand, the repulsive Coulomb interaction due to the presence of protons, 
even in a very small number, enlarges moderately the stability region of this matter in the 
range of 0.05 < p/po < 0.4, which is clearly shown in Fig. O Furthermore, Fig. [7] emphasizes 
and shows the important role of Coulomb interaction to stabilize matter without NT for 
almost zero perturbation. 

When the neutrino contribution is taken into account, the situation dramatically changes. 
This is shown in Fig. [8] for Yi e = 0.3 and in Fig. [9] for Y\ e = 0.4. For all parameter sets the 
instability boundaries expand, for which no substantial difference appears in the onset of 
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FIG. 9: (Color online) As in Fig. EJ but for Y h = 0.4. 

the instability due to the isovector adjustment, except for the region with q close to zero. In 
this region, the variation of the neutrino fraction in matter also yields an insignificant effect 
in the onset of instability, as we have expected. 

If we observe the longitudinal dielectric function at q close to zero (q = 0.01 fm -1 ), as 
shown in Fig. [TU1 then we can clearly see that the Z271* parameter set shows a quite different 
behavior compared to the other parameter sets, i.e., the Z271* parameter always yields El > 
in this limit. This fact indicates that the transition to more stable region at the points 
with small perturbation q < 0.1 fm" 1 around 0.1 < p/po < 0.3 of the Z271* is driven by a 
different mechanism compared to other parameter sets, i.e., it has a larger a sym that leads 
to larger proton and electron fractions (and less neutrinos for a fixed lepton fraction). As a 
consequence, a larger repulsion effect is produced compared to other parameter sets. 

From the fact that the adjusted parameter sets (indicated with *) have narrower instabil- 
ity region compared to their counterparts (indicated without *) for both models, then we can 
conclude that in this limit the onset of the instability is strongly related to the adjustment 
of the isovector sector. 

To dig up more information behind this fact, we show in Figs. [TT] and [12] the effects 
of the electrons absence (indicated by "— e" in the figure) and the absence of the Coulomb 
interaction (indicated by "—Coulomb) on the onset of the instability and on the instability 
in the limit of q close to zero, by switching off their contributions in the case of matter with 
NT. It can be seen from Fig. [TT] that the Coulomb interaction plays an important role in 
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FIG. 10: (Color online) Longitudinal dielectric function for matter with NT as a function of the 
ratio between baryon and nuclear saturation densities. All curves have been obtained by using 
the perturbation momentum q = 0.01 fm _1 and Y\ e = 0.4. The results are obtained by using the 
Horowitz-Piekarewicz (Z271 and Z271*) and Furnstahl-Serot-Tang (G2 and G2*) models. 
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FIG. 11: (Color online) As in Fig. El but calculated for the case of neutrino trapping with Yj e = 0.4. 



stabilizing the region with q < 0.6 fm _1 , in the range of 0.2 < p/po < 0.6. Furthermore, 
from Fig. [T5] we can see that if the electrons contribution were turned off, then the instability 
at q close to zero would disappear. Thus, for matter with NT the appearance of the large 
and negative El for the G2, G2* and Z271 parameter sets is caused by the fact that the 
repulsive interaction induced by the proton-neutron asymmetry (isovector) is unable to 
cancel the strong-attractive Coulomb interaction created by the presence of a substantially 
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FIG. 12: (Color online) As in Fig. [7J but calculated for the case of neutrino trapping with Yj e = 0.4. 

large number of electrons in matter. 

Figure [12] also shows that if the electrons are present but their Coulomb interactions 
were turned off (electrons behave as free particles) then the longitudinal dielectric function 
ei became larger and closer to El = but the instability boundary enlarges. This result can 
be used to emphasize the important role of the Coulomb interaction to stabilize the matter 
with NT in the limit of q close to zero. 

V. CONCLUSION 

We have studied how the instability region starts to appear in low-density matter de- 
scribed by the Horowitz-Piekarewicz and Furnstahl-Serot-Tang models. To this end we have 
utilized the longitudinal dielectric function at go = 0. The importance of the electron and 
Coulomb terms in matter with neutrino trapping has been investigated. It is found that the 
adjustment of the isovector terms has a more significant effect in the Horowitz-Piekarewicz 
model, i.e., producing a stronger repulsive isovector contribution which leads to a stronger 
correlation between its low density instability region and the a sym compared to the model 
of Furnstahl-Serot-Tang for matter without neutrino trapping. In the case of matter with 
neutrino trapping, the parameter sets with stiff EOS at low density lead to a large and 
positive £l(<1,<1o = 0). This demonstrates that, although the onsets of the instability of pa- 
rameter sets with stiff and soft EOS at low densities are similar, the driving mechanisms are 
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different. This fact might have an effect to the neutrino transport in matter. In both models 
the effect of the variation of the leptonic fraction is negligible, but the effect of the neutrino 
trapping on the onset of the instability region is significant. The presence of more protons 
and electrons in matter with neutrino trapping is the reason behind this phenomenon. The 
Coulomb term is found to be decisive in enlarging the stability of matter in this density 
region. The presence of the large and negative £z,(g, go = 0) in some parts of the instability 
region of matter with neutrino trapping originates from the fact that the isovector term 
is insufficient to cancel the attractive Coulomb interaction contributions generated by the 
presence of electrons (and protons). 
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